EXACT SOLUTION OF THE SIX- VERTEX MODEL WITH DOMAIN 
WALL BOUNDARY CONDITIONS. CRITICAL LINE BETWEEN 
FERROELECTRIC AND DISORDERED PHASES 



PAVEL BLEHER AND KARL LIECHTY 

Abstract. This is a continuation of the papers [1] of Bleher and Fokin and [5] of Bleher 
and Liechty, in which the large n asymptotics is obtained for the partition function Z„ of the 
six-vertex model with domain wall boundary conditions in the disordered and ferroelectric 
phases, respectively. In the present paper we obtain the large n asymptotics of Z„ on the 
critical line between these two phases. 



1. Introduction and formulation of the main result 

1.1. Definition of the model. The six-vertex model, or the model of two-dimensional ice, 
is stated on a square n x n lattice with arrows on edges. The arrows obey the rule that 
at every vertex there are two arrows pointing in and two arrows pointing out. Such rule is 
sometimes called the ice-rule. There are only six possible configurations of arrows at each 
vertex, hence the name of the model, see Fig. 1. 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



Figure 1. The six arrow configurations allowed at a vertex. 

We will consider the domain wall boundary conditions (DWBC), in which the arrows on 
the upper and lower boundaries point in the square, and the ones on the left and right 
boundaries point out. One possible configuration with DWBC on the 4x4 lattice is shown 
on Fig. 2. 
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Figure 2. An example of 4 x 4 configuration with DWBC. 



For eacli possible vertex state we assign a weight Wj, i = 1, . . . , 6, and define, as usual, 
the partition function, as a sum over all possible arrow configurations of the product of the 
vertex weights, 

Zn= j2 ^(^)' «^(^)= n^'^w=fi^^''^''^ (1-1) 

arrow configurations cr xdVn i=l 

where Vn is the n x n set of vertices, a{x) G {1, . . . , 6} is the vertex state of the configuration 
cr at the vertex x, according to Figure 1, and Ni{a) is the number of vertices of the vertex 
state i in the configuration a. The sum is taken over all possible configurations obeying the 
given boundary condition. The Gibbs measure is defined then as 

f.,.W = ^. (1.2) 

Our main goal is to obtain the large n asymptotics of the partition function Zn- 

In general, the six- vertex model has six parameters: the weights Wi. However, by using 
some conservation laws we can reduce these to only two parameters. Namely, first we reduce 
to the case 

Wi = W2 = a, ws = W4 = b, = Wq = c, (1.3) 
and then, by using the identity, 

Z„(a, a, 6, 6, c, c) = c^'z^ ( -, -, -, -, 1, 1 ) , (1.4) 

\c c c c J 

to the two parameters, ^ and ^ . For details on how we make this reduction, see, e.g., the 
works [1] of Allison and Reshetikhin, |12j of Ferrari and Spohn, or the work 



1.2. The phase diagram. Introduce the parameter 
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The phase diagram of the six- vertex model consists of the following three regions: the ferro- 
electric phase region, A > 1; the anti-ferroelectric phase region, A < — 1; and, the disordered 
phase region, — 1 < A < 1, see, e.g., [23]. In these three regions we parameterize the weights 
in the standard way: in the ferroelectric phase region, 

a = sinh(t — 7), 6 = sinh(t -|- 7), c = sinh(2|7|), < I7I < t, (1.6) 

in the anti-ferroelectric phase region, 

a = sinh(7 — t), 6 = sinh(7 + t), c = sinh(27), |t| < 7, (1.7) 

and in the disordered phase region, 

a = sin(7 — t), 6 = sin(7 + t), c = sin(27), |t| < 7. (1.8) 

The phase diagram of the six-vertex model is shown in Fig. 3. 




1 a/c 



Figure 3. The phase diagram of the model, where F, AF and D mark fer- 
roelectric, antiferroelectric, and disordered phase regions, respectively. The 
circular arc corresponds to the so-called "free fermion" line, when A = 0, 
and the three dots correspond to 1-, 2-, and 3-enumeration of alternating sign 
matrices. 



The phase diagram and the Bethe-Ansatz solution of the six-vertex model for periodic and 
anti-periodic boundary conditions are thoroughly discussed in the works of Lieb [20j-[23j, 
Lieb, Wu [21], Sutherland [27], Baxter [2J, Batchelor, Baxter, O'Rourke, Yung [3]. See also 
the work of Wu, Lin [29] , in which the Pfaffian solution for the six- vertex model with periodic 
boundary conditions is obtained on the free fermion line, A = 0. 

As concerns the six- vertex model with DWBC, it is noticed by Kuperberg [19j, that on 
the diagonal, 

^ = ^ = x, (1.9) 

c c 

the six-vertex model with DWBC is equivalent to the s-enumeration of alternating sign 
matrices (ASM), in which the weight of each such matrix is equal to s^-, where A_ is the 
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number of (— l)'s in the matrix and s = The exact solution for a finite n is known for 1-, 
2-, and 3-enumerations of ASMs, see the works by Kuperberg flQ] and Colomo and Pronko 
[9] for a solution based on the Izergin-Korepin formula. A fascinating story of the discovery 
of the ASM formula is presented in the book [7j of Bressoud. On the free fermion line, 7 = f , 
the partition function of the six-vertex model with DWBC has a very simple form: Zn = 1. 
For a nice short proof of this formula see the work [9] of Colomo and Pronko. 

Here we will discuss the ferroelectric and disordered phase regions, and we will use pa- 
rameterizations (11. 6 1) and (11.81) . The ferroelectric phase region consists of two connected 
components. For the sake of concreteness we will assume that 

7>0, (1.10) 

which corresponds to the component where 

b>a + c. (1.11) 

The parameter A in the ferroelectric phase region reduces to 

A = cosh(27). (1.12) 

1.3. Exact solution of Zn for finite n in ferroelectric phase. The six-vertex model 
with DWBC was first introduced by Korepin in [16j. In this paper, he derived an important 
recursion relation for Zn, which was subsequently used by Izergin [13] in deriving a determi- 
nantal formula for Zn in this model. A detailed proof of this formula and its generalizations 
are given in the paper of Izergin, Coker, and Korepin [13], see also the papers of Korepin, 
Zinn- Justin [18] and Kuperberg [T9] and the book of Bressoud [7]. When the weights are 
parameterized according to (II. 6p . the formula of Izergin-Korepin is 

_ [sinh(t - 7) sinh(t + 7)]"' 



where r„ is the Hankel determinant. 



+k-2. 



and 



, / N sinh(27) ,^ ^ ^, 

(^{t) = ^-TT tM^t r ■ 1-15 

^ ' smh(t + 7) smh(t - 7) ^ ' 

The aspect of the Izergin-Korepin determinantal formula that we exploit in this paper is 

that Tjv can be expressed in terms of related orthogonal polynomials, see the paper [30] of 

Zinn- Justin. In the ferroelectric phase the expression in terms of orthogonal polynomials 

can be obtained in the following manner. For the evaluation of the Hankel determinant, let 

us write 0(t) in the form of the Laplace transform of a discrete measure, 

ct^it) = . ""^^,^^.^1 r = 4^6-2*' sinh(27/). (1.16) 

^ ^ sinh(t + 7) sinh(t - 7) 

Then 

2 

J2 ^{k? n [2e-'*'" sinh(27/.)] , (1.17) 
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where 

l<i<j<n 

is the Vandermonde determinant. 

Introduce now discrete monic polynomials Pj (x) = + . . . orthogonal on the set N = 
{/ = 1, 2, . . .} with respect to the weight, 

w{l) = 2e-2" sinh(27/) = e-^^^^^' - e-^*'-^^', (1.19) 

so that 

oo 

J2Pjil)Pk{l)w{l) = hkSjk. (1.20) 

1=1 

Then it follows from (11.171) that 

n-l 
fc=0 

see Appendix in the end of the paper [6]. 



1.4. Exact solution of for finite n on the critical line between the ferroelectric 
and disordered phases. We consider the partition function Z„ on the critical line 

ha , ^ 
= 1. (1.22) 

c c 

We fix a point, 

aa— lha+1 

- = —.- = —; a.>l. (1,23) 
on this line, and we are interested in the large n asymptotics of the partition function 

/a — la — la + la + l \ 
Z.. = ^„(— . — . — . — .l.l). (1.24) 

Let us first derive a formula for Z„ on the critical line. To that end, consider the limit in 
the Izergin-Korepin formula in the ferroelectric phase, (11.131) . as 

t 

t,7^+0, >a. (1.25) 

7 



Observe that in this limit, 

a sinh(t — 7) a — 1 b sinh(t + 7) a + 1 

c ~ sinh(27) ' 2 ' c ~ sinh(27) ' 2 

By dH, 



;i.26) 



, a a 6 6 \ _ Znia, a,b,b, c, c) 

\ i-, i~ j — -2 ; U-^'j 



C C C C 



hence by f lTTT^ . and ffT^ . 



c c c c 



2 sinh()f: — 7) sinh(t + 7) 
sinh(27) 



-)2 n-l 



A:=0 ^ ^ 
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To deal with limit fll.25p we need to rescale the orthogonal polynomials Pk{l)- Introduce the 
rescaled variable, 

x = 2tl-2-fl, (1.29) 

and the rescaled limiting weight, 

a+1 



Waix) = hm (e 



-2tl+2'yl _ g-2ii-27^ 



a — 1 



> 1. 



;i.30) 



Consider monic orthogonal polynomials Pj{x; a) satisfying the orthogonality condition, 

Pj{x; a)Pk{x; a)wa{x)dx = hk,aSjk ■ (1-31) 

To find a relation between Pk{l) and Pk{x] a), introduce the monic polynomials 

Pk{x)=A''Pk{x/A), 



where 

A = 2t - 27, 

and rewrite orthogonality condition fll.20p in the form 

00 

J2Pj{lA)Pk{lA)wo,{lA)A = A^^+^hkS_ 



'jk, 



1=1 



which is a Riemann sum for the integral in orthogonality condition f ll.3ip . Therefore, 



lim Pfc(x) = Pk{x;a), 



and 



lim 



' ' '7 



Let us rewrite formula (11.280 as 

fa a b b 

\cccc 



2 sinh(t — 7) sinh(t + 7) 



sinh(27)A 

and take limit (ll.25p . In the limit we obtain that 



^2 n-l 

n 

fc=0 



(A;! 



\2 ' 



;i.32) 

1.33) 
;i.34) 

;i.35) 
;i.36) 

;i.37) 



a — 1 a — 1 a + l a + l 



1,1 



a + 1 



nhk,a 



fc=0 



1.38) 



Our main technical result in this paper will be the proof of the following asymptotics of 
hk^a- Let, as usual, 

C(s) = 1 + 7^ + tJ- + ..., Res>l. (1.39) 



Theorem 1.1. As k 00, 

hk,c 



In 



(A;! 



2^n{r - l)A;i/2 4/c 
A proof of this theorem will be given below 



^3/2^ 



a + 1 
a — 1 



1.40) 
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1.5. Main result. This paper is a continuation of the works [1] and [6], in which the large 
n asymptotics of Z„ is obtained in the disordered and ferroelectric phase, respectively. In 
[1] it is proven that in the disordered phase, for some e > 0, as n — oo. 



a a b b 



, , , ,1,1 =Cn^F- [l + Oin-% 
c c c c 



;i.41) 



where in parameterization (11.81) . 



27^ 



and 



12 37r(7r-7)' 

vr sin(7 — t) sin(7 + t) 
27 sin (27) cos 



nt 
27 



1.42) 



;i.43) 



The value of the constant C in (11.411) is not yet known. In [6] it is proven that in the 
ferroelectric phase, for any e > 0, as — > cxd. 



c c c c 



1 + e 



where in parameterization (11.61) . 



and 



sinh(t + 7) 
sinh(27) 



1.44) 



1.45) 



1.46) 



The main result of this paper is the following asymptotics of Z„ on the critical line between 
these two phases. 



Theorem 1.2. As n 00, 



a — 1 a — 1 a + 1 a + 1 



1, 1^ = Cn^G^F''' [1 + 0{n-^/^)] , (1.47) 



where C > 0, 



G = exp 



-c 



a — 1 

271 



and 



a + 1 



1.48) 



1.49) 
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The proof of Theorem 11.21 follows easily from Theorem II. 1[ Namely, from formula (11.381) 
and asymptotics (ll.40p we obtain that 



In 



7 ( ct—l a+l a+1 -| 1 \ 
I, 2 ' 2 ' 2 ' 2 ' ' / 



n-1 
fc=0 



n-l 



k=0 



h 



k,a 



1 



2A/7r(r - 1) A;i/2 4k 



-3/2N 



;i.5o) 



which implies Theorem 11.21 

1.6. Ground state configuration on the critical line. The ground state configuration 
a^'^ has the maximal weight. On the upper critical line between the ferroelectric and disor- 
dered phase regions we have that b = a + c. We also assume that a > 0, c > 0, hence a^^ 
should contain as many 6's as possible. The domain wall boundary conditions imply that in 

2 

each row there is at least one c. Therefore, any weight cannot be bigger than 6" ~"c". The 

2 

weight 6" -"c" does occur for the following unique configuration: 



; — 



Figure 4. A ground state configuration. 



if X is on the diagonal, 
(73 if X is above the diagonal, (1-51) 
(T4 if X is below the diagonal. 
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see Fig. 4, which is the unique ground state configuration with domain wall boundary con- 
ditions on the critical line. In the ferroelectric phase region, where 6 > a + c, the ground 
state configuration is obviously the same. If we set 



a — 1 , a + 1 



then the weight of the ground state configuration is equal to 

,2_ 



;i.52) 



w{a^') = F"^-", F = . (1.53) 



By (11.471) the ratio Zn/w{a^^) is evaluated as 

Z. 



■'n 



Observe that 



Cn'^G^F^fl + 0{n-^l% (1.54) 



limi^= limi^^^ = lnF, (1.55) 

SO that the free energy is determined by the free energy of the ground state configuration. 
This can be explained by the fact that low energy excited states are local perturbations of the 
ground state around the diagonal. Namely, it is impossible to create a new configuration by 
perturbing the ground state locally away of the diagonal: the conservation law N^i^a) = N4^{a) 
forbids such a configuration. Therefore, typical configurations of the six-vertex model in the 
ferroelectric phase region and on the critical line between the ferroelectric and disordered 
phase regions are frozen outside of a relatively small neighborhood of the diagonal. 

This behavior of typical configurations in the ferroelectric phase region and on the critical 
line between the ferroelectric and disordered phase regions is in a big contrast with the situa- 
tion in the disordered and anti-ferroelectric phase regions. Extensive rigorous, theoretical and 
numerical studies, see, e.g., the works of Cohn, Elkies, Propp [8], Eloranta [11], Syljuasen, 
Zvonarev ^26j, Allison, Reshetikhin pLj, Kenyon, Okounkov [15j, Kenyon, Okounkov, Sheffield 
[T7] . Sheffield [25], Ferrari, Spohn [12], Colomo, Pronko [10], Zinn- Justin [31], and references 
therein, show that in the disordered and anti-ferroelectric phase regions the "arctic circle" 
phenomenon persists, so that there are macroscopically big frozen and random domains in 
typical configurations, separated in the limit n ^ oo by an "arctic curve" . 



2. Large k asymptotics of h^^a 

We will use asymptotic formulae for orthogonal polynomials on [0,oo), obtained in the 
paper [28] of Vanlessen. To formulate and to apply the Vanlessen's asymptotic formula we 
will need to introduce some notations and to evaluate some parameters. Let us write 

Waiz) = e-' - e-'' = ze-'^^-'\ (2.1) 

so that 

Q{z) = z + \og^-^^, (2.2) 

where for log we take the principal branch with a cut at (— oo, 0]. Observe that the function 
Q{z) is analytic in a strip |Imz| < cq, cq > 0. Define the Mashkar-Rakhmanov-Saff numbers 
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A = as a solution to the equation 

— / Q'{^)\ ~s dx = k. 



27r Jo \j Pk-x 

As shown in [28], for large k there is a unique solution to this equation. 

2.1. Evaluation of f3k- By the change of variable, x = (3kU, equation fl2.3p reduces to 

^ Q'iPku)\ L ^ du= k. 



2vr Jo y 1 — u 

Set 

Oh = —r 
" 4k 



Then equation (12.41) reduces to 

2bk , , . I m" 

TT 

From (D, 



1 r - 1 

g'(^) = 1 + - 



z e('-i)^ - 1 ■ 
Observe that the function Q'{z) has poles at the points 

2rmii 

m = ±1,±2, 



r — 1 



After evaluating integrals of the first two terms of Q' (Ab^ku) , equation (12. 6p reads 



2k TT Jo e4('-"i)''fc'=" - 1 V 1 - M 
By the change of variable x = A;m, it reduces to 

^ ^ 2k nky^ io e^i-'^)''^- - 1 V 1 - (x/k) ^ 



Set 

1 



and consider the function, 



Observe that as e — 



l/e 



2 



r — 1)6 I X f°° {r — l)by/x dx 



2\ 



4(r-l)f>x _ 1 Y 1 _ ^2^, g4{r-l)fex _ ^ 



16^6(r - 1) 



( 



f(h^e) = b + / /: ^ J —dx-l, -<b<2. 



( 
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hence 

= (6 - 1) + ^ - ^ f^l^ + 0{e% (2.14) 
^ 8y 7ro(r — 1) 

It is easy to see that equation fl2.13p can be differentiated in h infinitely many times, and 
hence the function f{b,e) is C°° in a neighborhood of the point 6=1, e = 0. In addition, 

/(1,0) = 0, ^^^ = 1- (2.15) 

By the implicit function theorem, there is a C°°-solution b{e) of the equation f{b,6) = 0. 
From fl2.14p we obtain that 

6(e) = l-- + -^4M= + 0(.^), 5-0. (2.16) 
2 8^7T{r - 1) 

Since = b{k~^^'^), this gives 
By (ESD, 

Pk = 4:kbk. (2.18) 

2.2. Evaluation of the equilibrium measure. Set 

Vkix) = ^Q{Pkx), (2.19) 

and consider the following minimization problem: 

E = inf J(/i), (2.20) 

where 

/(yu) = - jj log |x - y\ dfi{x)dfi{y) + y Vk{x)dfi{x). (2.21) 

and inf is taken over all probability measures on [0, oo). There exists a unique minimizer, 

/i = /ifc, and it has the following properties: 

(1) The support of /i^ is the interval [0, 1]. 

(2) The measure /i^ is absolutely continuous with respect to the Lebesgue measure. 

(3) The density function of /i^ has the form. 



where qk{x) is analytic and positive on [0, 1]. 

The equilibrium measure /ifc is characterized by the Euler-Lagrange variational conditions: 
there exists G M such that 



2 / log \x - y\dfik{y) - Vk{x) -lk = for x G [0, 1], 
Jo 

2 [ log la; - y\diikiy) - V^ix) - 4 < for a; ^ [0, 1]. 

^0 



(2.23) 
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The function qk{z) in (12.221) is given by the formula, 

where \J is taken on the principal branch, with cut on [0, 1], and F is a positively oriented 

contour containing [0, 1] U {z} in its interior, with the additional condition, that the function 
Vl{y) is analytic inside T. By flCT]) and Q, 

V,{z) = Ahz + - In (2.25) 



hence 

where 
hence 



Qkiz) 



+ (2.26) 



7fc = 4(r - l)bk, (2.27) 

dy 



2711 /p y ?/ — 1 



46. + ^- 



ky e^kky _ 1 



y- z 



By taking the residue at infinity, we obtain that 



(2.28) 



q,iz) = Ah + s,iz), s,{z) = -^.i ,/^— T-^^^T r • (2-29) 



2vri /r Y y - 1 {e^'^^y - l){y - z 
Observe that the function Vl{z) has poles at the points 

2 77" 772 

Zn = ^, r2 = ±l,±2,..., (2.30) 

hence the contour F has to pass close to 0. We choose F such that 

^ > dist(0, F) > ^, ci > dist(l, F) > C2 > 0, (2.31) 
k k 

see Fig. 5. More precisely, let for a given 2; G C, m{z) G [0, 1] be the closest point from z on 
[0, 1], so that 

inf{|z - m|, M G [0, 1]} =\z- m{z)\. (2.32) 




Figure 5. The contour F. 
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Then we define, for a given 6 > 0, 



r = T{5, k) = {zeC: \z- m{z)\ = 5 



1 

- + m[z) 



(2.33) 



and we clioose 6 to be sufficiently small so that the points Zn in (12.301) he outside of F, see 
Fig. 5. Observe that F(0, k) = [0, 1]. 

With the help of the change of variables, u = ky, we obtain that 



^'^''^ 2inky^ ip J (u/k) - 1 (e^^« - 1)(m - kz) ' ^^'^^^ 



which implies that 



or even that 



sup \sk{z)\=0{k-^/^), (2.35) 

0<z<l 



sup sup \skiz)\=0{k~^/^). (2.36) 

o<d<^ zer{d,k) 



For z > 1, the function Sk{z) can be reduced to 



''^'^ ~ ^'^1 z- 1 (e^^^^ - 1) ? / {ef^'^y - l){z - y) ' ^^'^'^^ 



It implies that 



where 



Sk{z) = — + rk{z), (2.38) 



-^ = --1 x[^^4^. (2-39) 



and rfc(z) satisfies the estimate. 



\r,{z)\ < ^£^— z>l, (2.40) 
z\/ z — 1 fc^/"' 



with some C > 0. Indeed, 



/ ^ 1 Ik ( I y ydy . . 

'■'^''^ ^'V -2 - 1 e^^'^ - 1 vr io V 1 - 1/ (e^^'^^ - 1)(^ - 1/)^ ' 

The first term on the right is exponentially small in z and k, and it obviously satisfies 
estimate fl2.40p . In the second term on the right, let us split the integral in two integrals, 
from to I and from | to 1. The first part is estimated as follows: 



(e7.fcy_l)(^_y)^ - ^2;^ ^e^,ky_i^ 
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hence it satisfies estimate (12.40p . For the second part we have that 



< 



y 



ydy 



< 



1 r 

V2 _ 1) J 



dy 



1 - y {e-y^^y -l){z- y)z " z(eT'=*^/2 - 1) ji (z - y)^/T^ 

1 du ^ Ci 



(2.43) 



< 



u = l-y, 



which satisfies estimate (12 .40^ . Thus, (12.401) is proved. 
From (I2.38P we obtain that 

gfc(z) =46fc + ^ + rfc(z), (2.44) 

where is given by formula (12.391) and rk{z) satisfies estimate (12.401) . 
2.3. Evaluation of the Lagrange multiplier. Introduce the function 

gk{z) = / \og{z - x)dij,k{x), zeC\[0,l], 
Jo 

where the branch of log is taken on the principal sheet, with a cut on {—oo, 0]. Then 



Mz) = g'kiz) 



z — X 



z e C \ [0, 1]. 



From equation (12.451) it follows that as 2; — oo, 

Qkiz) = log z + 0{z-'] 

and from (I2.46p . that 



UJk{z) = - + 0{z' 



From equation (12.230 it follows that 



UkiZ 



Vliz) z-lqk{z) 



2 y z 2 
see, e.g., equations (3.27), (3.29) in [28], and 

h = 2gkil)-Vkil). 

Since 



nu 

-gk{l) = lim [gk{u) - log u - gk{l)] = lim / 

U—fOO U—rCO 



UJk{z) 



dz 



Vl{z) lz-l_ qkjz) 1 
z 2 z 



(2.45) 

(2.46) 

(2.47) 
(2.48) 

(2.49) 
(2.50) 

(2.51) 



dz. 



we obtain that 



hence 



^^^■(1) 



V^iz) I Z -I Qkjz) 1 

z 2 z 



dz. 



z-l 2 

QkKz) - - 

z z 



dz-Vk{l). 



(2.52) 



(2.53) 



EXACT SOLUTION OF THE SIX- VERTEX MODEL 

We split the last integral as 

' 2 



z- 1 



dz 



qk{z) H 46fc - — 

Z Z KZ 



V;'(z)-46,-— dz 



dz = Ii — I2. 



From fl2.26p we have that 



h 



Vl(z) - Ah 



kz 



dz 



Ikdz 
1 eTfe''^ - 1 



0{e-'''^) , Co > 0. 



Let us evaluate h- By f l2l44|) 

h-- 



-2 — 1., Z — ^ CLk Z — 1 2 1 

Abk + d - + \ rk{z) + - ~ Ah ~ — 

z \ z z \ z z kz 



dz. 



Since 



z 2z 



dz = In 2, 

2 



we obtain that 

/2 = 6fc(2-41n2) + 



2; - 1 ttfc z - 1 , . 2bk-2 1 

+ \ rk{z) — 

z z \ z z kz 



dz. 



From estimate (12.401) we obtain that 



z-1 



n{z)dz = 0(r^/2), 



hence 



/2 = 6fe(2-41n2) 



z-1 at 2bk-2 + l 



z z 



dz + 0{k-^^^). 



From equation (12.491) we have that 



Abk + — + rkiz) 
z 



By equating terms of the order of - on both sides, we obtain that 



hence 



' = 2^-y + ^^' 



Ofc = 26fc - 2 + ^ . 



By substituting this expression into (I2.60p we obtain that 

.00/ ^ 



J2 = 6fc(2-41n2)+ 26fc-2 + 



k 



z 



1 \^ + o{k-'/'). 

z 
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Since 



we obtain that 



z - 1 



dz 



21n2 - 2, 



(2.65) 



/2 = 6fc(2-41n2)+ ( 26fc-2 + ^ ) (2 ln2 - 2) + 0(A;~^/2^ 



21n2-2 

-26fc - 4 In 2 + 4 + ^ + 0{k~^l^) 



k 



(2.66) 



= 2-41n2 

By (ESI, (EMD, 

and by fl2:25l) and fl2T7D . 



21n2 - 1 

k 



/fc = /2-/l-V^(l), 



+ o(r^/2). 



(2.67) 



W)=46. + i^ + 0(rV^) 



, InA; 21n2-2 C(|) 1 

4 H ^ ^ 

k k 2^7f{r - l)A;3/2 2fc^ 



+ 0(r^/2) 



(2.68) 



hence 



Ink 1 
-2-41n2 — + T 



3C( 



+ TT77 + 0(^-'/') 



A; k 4^7T{r - l)A;3/2 2^2 
2.4. Evaluation of /ifc^Q. According to Vanlessen's asymptotic formula, see 



(2.69) 



^ n2k+2kU 



1 + 



3 47 

+ 



4g,(0) 12gfc(l) 4gfc(l)2; A; 



- + 0(fc- 



(2.70) 



Observe that in [28j| this formula is proved under the assumption that the weight for orthog- 
onal polynomials has the form 



w(x) = x'^e "^^""^ 



7 > -1, 



(2.71) 



where Q{x) is a polynomial. In Appendix at the end of the paper we show what changes in 
the paper of Vanlessen [2^ should be made to prove f l2.70p for the weight Wa{x) given by 
formula By 



By draD, 



hence 



gfc(O) = 4 + 0(fc-V2)^ ^^(1) = 4 + o(fc-V2). 



7fc 



27riA;i/2 J^^ y (^^/A;) - 1 {e<ku -!)(«- A;)2 ' 



(2.72) 

(2.73) 
(2.74) 
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because by condition (12.311) . the function is bounded by 1/c for m G kV. From (12.721) 
and (I2.74P we obtain that 

By substituting formulae fl218D . (12X71) . fl239D . and fl2T5|) into fl2T0|) and by using the 
Stirhng formula for A;!, we obtain that 

In ^ = ^ilL + -L + 0(A;-3/2) ^2.76) 

(we use MAPLE for this calculation). Theorem 11.11 is proved. 

Appendix A. Proof of formula (12.701) 

We use the notations and results from the work [28j of Vanlessen. The essential difference 
with [28] is that we consider not a fixed but a shrinking neighborhood of the origin, 

U6,k=^zeC: kl<^}, (A.l) 

where 5 > is small enough so that the function Vfc(x) is analytic in t/^^, see (I2.19p . 
As in [28j, we consider a sequence of transformations of the Riemann-Hilbert problem for 
orthogonal polynomials, and in the end we arrive at the following Riemann-Hilbert problem 
on a 2 X 2 matrix- valued function R{z): 




Figure 6. The contour Tr. 

(1) R{z) is analytic on C\r/{, where is the contour shown on Fig. 6, and it has limits, 
R+{z) and R-{z), on P/j, as z approaches a point on Vr from the left and from the 
right of the contour, with respect to the orientation indicated on Fig. 6. 

(2) On P/j, R{z) satisfies the jump condition, R+{z) = R^{z)vr{z), where vr{z) is an 
explicit matrix-valued function. 

(3) R{z) ^/+^ + ...asz— > cxD, where / is the identity matrix. 

The contour Vr consists of the circle dUs^k, the circle dUs, where 

Us = {zeC: \z-l\<6}, (A.2) 
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the boundaries of the lenses, 7-1- , and the semi-infinite interval [1 + S, 00). The jump matrix 
vr on dUs^k has the following asymptotics as — 00: 



vr{z) 



(A.3) 



n=l 



see formulae (3.105) and (3.98) in [28]. This asymptotic formula holds under the condition 
that k'^z —>■ 00. Under this condition, for any > 1 there exists a constant Cn > such 
that 

N 

Vr{z) - I - Y,~^n{z)k-^ 
n=l 

The condition k'^z ^ cxd is valid for z G dUs^k, and in this case the last estimate gives that 



< 



C 



N 



(A.4) 



sup 



N 



vr{z) - I - ^ A„(z)A;" 



n=l 



< 



c 



N 



k 2 



C 



C 



N 



N 



The coefficients A„(z) in (1A.3P are given by the following formula: 



K{Z) = ^^p(-)(^)(-^)|-3^J_^)-|.3p(oo)(^)-l^ 
0fc(2;)«/2 



where 



a{z)+a(z) a{z)—a(z) ^ 

2 ~" 2i 

a{z)—a(z) a{z)+a{z) 

2 



'2z-l + 2^z{z - 1) 



^-3 



a{z) 



z - 1 



1/4 



0fe(2;) = - 



1 - S 



-\ 2 



qk[.s) ds 



n;^i[4-(2j-l)V (_ll::(3 + 2 



16"ri! 



-l)"+i(n-i)z i-(3 + 2n) 



(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 



see equation (3.99) in [28j. The function 4>k{z) is analytic in Us^k- From fl2.29p . fl2.17p . and 
f l2.36p we obtain that as /c — > 00, 

(A.IO) 



sup \qk{z)-4\ = 0{k-'/^). 



By (1A.8P this implies that 



sup 



<Pk{z) 



0{k 



-1/2N 



1 - S 



4^2 

ds\ =-4z + — + ... (A.ll) 



The function A„(z) is meromorphic in f/^^fc with the only possible pole at the origin of the 
order at most [^^] , see [28]. This result, combined with explicit formula (lA.6p . implies that 



there exists c„ > such that 



sup 

zedUs^k 



K{z)k- 



-n+ 



n + l] 



(A.12) 
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This, in turn, allows us to improve estimate (1A.5I) as follows: for any > 1 there exists 
Cat > such that 

N 

vr{z) - I - ^ A„(2;)A;" 



sup 



n=l 



(A.13) 



so that the error term is estimated by a constant times the estimate of Ajsf+i{z)k 
When N = 1, this gives that 



sup 



vr[z) 



k 



Oik-'] 



The function Ai(2;) has a simple pole at and its residue is equal to 



Bi, 



-(73 



1 i 
i -1 



(A. 14) 



(A.15) 



16g,(0) 

see equation (4.11) in [28]. The function Ai(z) — — is regular at z = and from explicit 
formula (]A.6I) we obtain that as /c — *• oo, 



sup 

z(^dUs.k 



Z 



0(1] 



hence from flA.14p we obtain that 



sup 



Vii{z) - I -y- 
kz 



0{k 



(A.16) 



(A.17) 



The problem here is that vr{z) is not close to I on dUs^k- We will overcome this obstacle by 
a transformation of the Riemann-Hilbert problem for R{z). 
Observe that 

Tr5fc = 0, detSfc = 0, (A. 18) 

hence 

= 1, (A.19) 



det(/ + ^ 

kz 



hence the matrix / + is invertible for any z ^ 0. Let us make the substitution. 



R{z) 



R{z), z e fJs,, 



(A.20) 



Then R{z) solves the Riemann-Hilbert problem similar to the one for R{z), with the jump 
matrix vpi{z] such that 

-1 

z e dUs,k (A.21) 



Vr{z) = Vr{z) { I 



and 



vr(z 



Bu 



kz 



I + t^\vr{z)\I+t^ 



B, 



kz 



(A.22) 
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From ( ]A.17[) we obtain that 



sup \vji{z) — /| = 0{k 



(A.23) 



Also, since the equihbrium density function diverges as z at the origin, we obtain that 
vr{z) is sub-exponentially small on the boundary of lenses. 



0(e-"^), 

This implies that vr{z) satisfies a similar estimate, 

sup \vR{z)-I\=0{e-'^), 

In addition, 



sup \vr{z) - /| 
2G74. U7- 



0. 



0. 



sup \Vr[Z) 

z&aus 



and 



dniz)- l\=0{e-'"''), z>l; c> 0. 



(A.24) 



(A.25) 



(A.26) 



(A.27) 



These estimates of smallness of {vji{z) — J) on Tr enable us to solve the Riemann-Hilbert 
problem for R{z) by a series of perturbation theory. The fact that the radius of Us,k^ r — ^ 



is tending to zero does not cause a problem, see appendix to the work [5J of Bleher and 
Kuijlaars. 

The rest of the proof of formula fl2.70p goes along the lines of [28]. Namely, by formula 
(4.17) in 



hk,a = I (3f+^e^'' [I - m{Ri)i2 + 0{k~ 



By flCT]) . 



(i?l)l2 = (i?l)l2 + + 0{k-^) = + 



3i 



By applying formula (4.11) in [28j to -Ri, we obtain that 



+ 0{k-^) 



(A.28) 



(A.29) 



I 12 



+ 



47i 



QAqkilfk 192g,(l)A; 



+ 0{k- 



(A.30) 



Observe that the first term in formula (4.11) in [28j is missing in this case, because the 



function 



A 1(2;) — ^ is regular at 2; = 0. From the last two formulae we obtain that 



1 12 



^ 47 



_4gfc(0) 4gfe(l)2 12^^(1) 
By substituting this into flA.281) we obtain fl2.70p . 



^ + 0(0 



(A.31) 
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